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Abstract 

The Random Matrix Model approach to Quantum Chromodynamics (QCD) with non- vanishing 
chemical potential is reviewed. The general concept using global symmetries is introduced, as well 
as its relation to field theory the so-called epsilon regime of chiral Perturbation Theory (exPT). 
Two types of Matrix Model results are distinguished: phenomenological applications leading to 
phase diagrams, and an exact limit of the QCD Dirac operator spectrum matching with e%PT. All 
known analytic results for the spectrum of complex and symplectic Matrix Models with chemical 
potential are summarised for the symmetry classes of ordinary and adjoint QCD, respectively. 
These include correlation functions of Dirac operator eigenvalues in the complex plane for real 
chemical potential, and in the real plane for imaginary isospin chemical potential. Comparisons of 
these predictions to recent Lattice simulations are also discussed. 
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1 Introduction 



In the phase where chiral symmetry is broken and quarks and gluons are confined Quantum Chromo- 
dynamics (QCD) is strongly coupled and perturbation theory breaks down. Several methods alow to 
study this situation from first principles. In an effective field theory approach only the lowest degrees 
of freedom are taken into account, such as the pseudo Goldstone bosons in chiral Pertubation Theory 
(xPT). Second, in Lattice Gauge Theory (LGT) QCD is studied on a Euclidean space-time lattice. 
At zero baryon chemical potential [i = LGT has been very successful, e.g. in describing the particle 
spectrum or the chiral phase transition. Turning on \x in the phase diagram, LGT faces the so-called 
sign problem: the action becomes complex in general. This invalidates the usual method to weight 
configurations by their Boltzmann factor, as it ceases to be a real positive quantity. 

In this article we will study a simple model called chiral Random Matrix Model (MM) in order 
to understand some aspects of this problem. This includes insights into chiral symmetry breaking at 
non-zero chemical potential, and quantitative aspects of the Dirac operator spectrum close to zero 
momentum, illuminating the interplay between LGT and xPT. MMs can be mapped to field theory 
in a particular limit called the epsilon regime or exPT. Clearly QCD is not only a MM, and we will 
cover aspects dictated by global symmetries alone. Also our MM does not solve the sign problem on 
the Lattice, but it provides a model that can be exactly solved despite this problem, and that offers 
quantitative predictions for the Lattice. The three types of real, complex and symplectic chiral MMs 
match with the classification of gauge theories in general (including QCD) into 3 classes, displaying 
different ways or the absence of the sign problem. 

To start let us briefly summarise what has been achieved for MMs of QCD at zero chemical potential 
since their introduction in pQ. An excellent review describing its concepts and success can be found in 
[2]. The introduction of MMs as a tool to study complicated physical systems by means of symmetry 
predates QCD substantially, going back to Dyson and Wigner in the late 1950 and early 1960. Initially 
developed for heavy nuclei MMs have today found applications in most fields of physics (and other 
sciences), as reviewed in [3]. The first MM for the Dirac spectrum of QCD was written down and 
solved in pQ, including all quenched and massless eigenvalue correlation functions for the MM Dirac 
operator eigenvalues. It became apparent that the 3 fold classification of Dirac operators for different 
field theories according to charge and complex conjugation symmetry [5j exactly coincides with the 3 
possible chiral symmetry breaking patterns [5]. The corresponding 3 MMs are called chiral ensembles, 
having real, complex and quaternion real matrix elemental]. All eigenvalue correlation functions for 
these 3 chiral MMs ensembles were consecutively computed, including an arbitrary number of massive 
quark flavours N f [H El [9j Q3H 031 021 Q3] , and individual eigenvalue distributions [Til HP] . 

The equivalence of such MMs to e^PT was first realised on the level of partition functions [Tj I15j. 
leading to identical sum-rules for the Dirac operator |16t 117]. As explained by Gasser and Leutwyler 
|18] the e^PT regime also called microscopic domain is reached, when the zero-momentum Goldstone 
modes completely dominate the partition function, decoupling from Gaussian propagating modes. 

The MM density [19] and density-density [20] correlations were shown to coincide with those 
derived from e^PT alone, and consequently the individual eigenvalue distributions to leading order 
[21] . Note that so far no higher density correlations, nor density correlations in symmetry classes other 
than QCD have been computed from e^PT. 

Furthermore it was understood precisely when the MM-e^PT correspondence breaks down, at the 
so-called Thouless energy [22], making the MM approximation transparent. Even when going beyond 
the MM regime its techniques have been useful to compute expectation values of vector and axial 
vector currents, see [23] and references therein. 

1 Non-chiral MMs as those of Dyson and Wigner can be used to study flavour symmetry breaking in 3 dimensions, 
see [6] [7] and references therein. 
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Below the Thouless energy MMs give precise analytic predictions for the low Dirac operator eigen- 
values lying inside the e-regime, as a function of quark masses and gauge field topology, for each of the 
3 symmetry classes. Such quantitative predictions have been extensively checked in all 3 classes, in 
quenched and unquenched simulations. First tests were done using staggered fermions, giving confir- 
mation in all 3 symmetry classes [24] in the sector of topological charge v = 0. For staggered fermions 
care has to be taken compared to the global symmetry in the continuum in 2 of the 3 classes. 

The predictive power of MMs has gained even more relevance in the important development of 
chiral fermions, that possess an exact Lattice chiral symmetry and well defined topology away from 
the continuum (see [25J for a review). Here MMs have been used to test such algorithms at fixed 
gauge field topology, being confirmed for all 3 classes by various groups [26] . Subsequently staggered 
simulations have been improved, also showing the correct topology dependence as overlap fermions 
[27]. 

In this article we will try to go through the same steps for MMs with non-zero chemical potential, 
and compare the analytical prediction available for 2 of the 3 classes to Lattice data whenever possible. 
Very recently a review on the same topic appeared [28], focusing more on the relation to e^PT than 
the MM concepts, its solution and comparison to LGT covered here. 

In order to emphasise the limitations of our approach let us recall some important questions of 
LGT at finite density and temperature T. One of the purposes underlying Lattice simulations is to 
better understand the restoration of chiral symmetry, going hand in hand with deconfinement for 
QCD. Major experiments such as RHIC as well as the transition in the early universe happen at very 
small fj, close to the temperature axis. Here several techniques have been applied to avoid the sign or 
complex action problem, where we refer to [29] for a review and [30] for recent developments. Physical 
questions as for the free energy, pressure, equation of state, screening mass or string tension have been 
addressed, and the dependence on the number of flavours and quark masses is being tested (see [31] 
for a review and references). However, none of these techniques actually solves the sign problem. For 
that reason the phase diagram far away from the T-axis is still current subject of debate [32] . 

Our simple MM cannot give answers to these physical questions, despite providing a qualitative 
(flavour independent) phase diagram for QCD. The domain of our MM solution is at small or zero T 
and small to moderate \x. Its analytical predictions there can serve as a testing ground, both in the 
presence or absence of the sign problem. It also can tell us in which region simulations are doable 
at finite volume and when the sign problem hits in, as reviewed in [28]. Because of this situation so 
far we can only provide evidence for MM predictions by comparing to quenched Lattice QCD, and to 
unquenched LGTs not suffering from the sign problem. 

Several topics using MMs at finite [i are not covered in this review. In section [3] we review 
only aspects of the phase diagram for QCD. Similar results exist for two colours [33] as well, and 
for a discussion of the chiral Lagrangian with chemical potential for two colours and the adjoint 
representation we refer to [33]. The phase diagram of QCD can also be analysed using the complex 
zeros of the partition function, where we refer to recent [35] as well as previous work [36]. Recently the 
phase of the fermion determinant has been computed [37], which can be tested in Lattice simulations. 
In addition to chemical potential and temperature also colour degrees of freedom can be taken into 
account in MMs [38] in order to model single-gluon exchanges, or the symmetries of the staggered 
operator can be added [39]. However, these additions typically destroy the exact analytical solvability. 
Last but not least a factorisation method has been advocated by the authors of [40]. It uses MMs 
as a toy model in order to tackle the sign problem in Lattice simulations. As mentioned earlier MMs 
enjoy many different applications in Physics. For applications of MMs with complex eigenvalues in 
other areas we refer to the review |41| . 

This article is organised as follows. In the next section [2] we explain the concept of our MM 
approach including a symmetry classification and its relation to exPT. After a brief view on the 
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resulting phase diagrams in section [3j we turn to the relation between Dirac operator eigenvalues 
and chiral symmetry breaking in section [U Here we summarise what is know for complex eigenvalue 
correlation functions in the class containing QCD, and confront with quenched lattice results. Section 
[5] contains all results for symplectic MMs corresponding to field theories with fermions in the adjoint 
representation. Here we compare to unquenched Lattice results. In section [6] we present the solution 
of a MM for the QCD class with imaginary isospin chemical potential where unquenched simulations 
are straightforward. Conclusions and an outlook are given in section [71 

2 Concept of Matrix Models at Finite fi 
2.1 Global symmetries 

The general strategy in applying Matrix Models to Physics is to replace a complicated operator by a 
random matrix with the same global symmetries, trying to gain some information on its spectrum. In 
1993 Shuryak and Verbaarschot applied this idea for the first time to the Dirac operator Jj) in QCD 
[I], which was then extended to QCD with chemical potential by Stephanov [42J: 

*M = * + _,,%„ < 21 > 

The matrix $ has constant, that is space-time independent matrix elements, typically with identical 
Gaussian distributions and hence its name random. Here we have used already 2 global symmetries: in 
the continuum chiral symmetry implies that {Jp([x), 75} = 0, leading to an off-diagonal block structure. 
Second, we have used that the Euclidean Dirac operator is anti-hermitian at zero chemical potential, 
and that /ijq is hermitian. 

The average over the gauge fields of the Yang-Mills action is replaced by a random matrix average, 

expl-Sgauge] — > exp[-a 2 Tr $t$] , ( 2 . 2 ) 

when computing observables. The inverse variance a 2 in the Matrix Model is to be fixed later. Using 
a lattice regularisation in a finite volume V the Dirac matrix becomes finite dimensional, and we thus 
can choose the random matrix of dimension N x [N + u). Here we have furthermore restricted 
ourselves to a sector of fixed topology u, or through the index theorem to a fixed number of zero 
modes of Ip. Its eigenvalues are defined as 

= iz k (p k , k = l,...,V (jj, = => z k E M) . (2.3) 

For zero fi they lie on the imaginary axis, whereas for nonzero fi the operator lp(fJ,) has complex 
eigenvalues in general, being a complex non-Hermitian operator. The large- iV limit in the matrix size 
thus corresponds to the thermodynamic limit lim2iV ~ V — > 00, keeping the size of the gauge group 
Nc fixed. 

What are the matrix elements of <£? It turns out that they can be further restricted by symmetry 
and fall into 3 different classes. According to the choice of gauge group and its representation it 
was shown for gauge groups SU(Nc) at fj, = [4J and at fj, ^ [M| that for either two colours 
Nc = 2 or the adjoint representation an anti-unitary symmetry allows to choose the matrix elements 
&ij to be real or quaternion rea@, respectively. For the class with SU (Nc > 3) in the fundamental 
representation containing QCD no such symmetry exists, and the matrix elements are complex. These 
are the global symmetries that the random matrices and the exact Dirac operators share, as shown in 
the table [TJ This classification can of course be extended to any gauge groups having a real, complex 
or pseudo-real representation^, and we have added U(l) to the table. 

2 For quaternions 1| denotes the quaternion unity matrix. 
3 I would like to thank P. Damgaard for this comment. 
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SU{2) fund. 


5i7(iV c > 3), 17(1) fund. 


SU(N C ) adjoint 


Th( \ ( $ + ^\ 
^={-* T + ^ ) 


/ o $ + % A 

V -& + 1\fJL /" 


/ $ + % \ 
V — *t + ij M o y 1 


Qije R,p D = l 


* -€ C,P D = 2 


real HI , (3 D = 4 


as [Ca 2 K, 0] = 
3 finite fraction zj~ G R 




as [CivT, 7/>] = 
Mk: z k ^ R, iR 


SU(Nc) adjoint staggered 


SU{N C > 3), 17(1) fund, 
staggered 


SU (2) fund, staggered 



Table 1: Table of global symmetry classes. 



The anti-unitary symmetry in the real and pseudo-real case is provided by a product of the charge 
C and complex conjugation matrix K, and a Pauli matrix 02 for SU(2) [UGH]. We have introduced 
the Dyson index (3d to label the number of independent degrees of freedom per matrix element of <£. 
In the class (3d = 1 the non-symmetric, real matrix Ip(fi 7^ 0) has a real secular equation. Hence it can 
have real eigenvalues^] in addition to complex eigenvalues occuring in conjugated pairs. For that reason 
this class has a real det [$(//)] which is not necessarily positive. This sign problem is thus reduced 
from a phase to a fluctuating sign. For an even number of two-fold degenerate flavours the product of 
all determinants remains positive and has no sign problem. For (3d = 4 the determinant det [$(//)] is 
alway positive, as the eigenvalues of a real quaternion matrix occur in complex conjugated pairs [44) . 

Notice that for discretised Dirac operators on a Lattice their global symmetries may change: when 
using staggered fermions the symmetry classes (3d = 1 and 4 have to be interchanged [15]. On the 
other hand Wilson type Dirac operators are in the same symmetry class as the continuum. So far 
the transition from the discrete to the continuum symmetry class has not been observed for staggered 
fermions in the Dirac spectrum. 





/x = m = 




(3d = 1 


SU{2N f ) -» Sp{2N f ) 


SU(N f ) L 


x SU(N f ) R -» 


Sp(JV/) L x Sp(JV» fl 
Sp(Nf)v for m / 


(3d = 2 


SU{N f ) L x SU(N f ) R -» SU{N f ) v 


SU(N f ) L 


x SU(N f ) R -» 


5Z7(iV / ) v 


(3d =4 


SU{2N f ) - SO(2N f ) 


SU(N f ) L 


x SU(N f ) R - 


SO(JV» £ x 50(iV/) fl 
SO(N f ) v for m / 



Table 2: Table of spontaneous (left) and explicit (right) chiral symmetry breaking patterns. 



4 In a different MM their expectation vaiue growns ~ viV [43] . 
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It turns out that the classification into 3 groups exactly corresponds the 3 possible patterns of 
spontaneous chiral symmetry breaking introduced by Peskin [5], where a maximal breaking is assumed. 
Here we have dropped both the vector UyiX) and axial vector factor as they are not affected 

by the spontaneous breaking. The case of symmetry breaking in 2 dimensions is special, and we refer 
to [36] for an exhaustive treatment of the Schwinger model and its relation to MMs. 

The mass term in ( see ec l- (|2.1U|) ) of course explicitly breaks chiral symmetry. At \i = 

the patterns are the same as for spontaneous breaking in all 3 classes, see table [2j A nonzero /U-term 
also explicitly breaks chiral symmetry, and the breaking patters differs from the spontaneous one for 
(3d = 1,4 as discussed in [34j . both for m = and m ^ 0. There is a competition for the vacuum 
aligning to the chemical potential or the mass term if m ^ 0, and we refer to [34] for details. Only in 
the (3d = 2 class the explicit and spontaneous breaking patters remain the same at fj, ^ 0. 

This complication for (3d = 1,4 at \x ^ is to be contrasted with the global symmetries in table 
[TJ Despite breaking the anti-hermiticity /i / does not change these global symmetries. Therefore 
in the following we use tabled] as a guiding principle to build our MMs with /i ^ 0, following [42] for 
(3d = 2 and [47J for (3d = 1,4 where these MMs were first introduced. A second reason is that the 
effective field theories for (3d = 1 and 4 are much complicated due to the structure of the coset space 
for the Goldstone manifold from table [2] Only limited analytical results are available even at \i = 0, 
see pH US]. 

]p(n) real JPip) complex Ip{p) quaternion real 
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Figure 1: Scatter plots from [17] of complex eigenvalues in the 3 different MM symmetry classes 
(columns) obtained numerically for a given JV>1 and two different values of \i. The larger \x (lower 
row) is just below the critical fx c ~ 0.527 in these models. 

In [47] the effect of chemical potential in the three different symmetry classes in table[]]was analysed 
using a MM, as shown in fig. [TJ It shows a scatter plot of complex MM eigenvalues obtained numerically 
for large random matrices. In the matrix model (as well as on the Lattice) the eigenvalues stay inside 
a bounded domain. The shape of this boundary is a function of fi and was computed in a saddle point 
calculation at N — > oo for a non-chiral MM [42] for a 2 = N: (y 2 + 4/i 2 )(^i 2 - x 2 ) 2 + x 2 = 4/x 2 (/i 2 - x 2 ) 
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where z = x + iy. This envelope is plotted as a full line in fig. [TJ It splits into two sets at [i = 1 a 
spurious transition which is beyond the chiral transition at fi c ~ 0.527 for T = 0, see eq. ()3.3p . The 
envelope generalises the Wigner semi-circle for the MM spectral density at /i = 0, which may also 
split into two arcs due to the influence of an external field like temperature T [49] . 

Fig. Q] illustrates the signature of the 3 symmetry classes: for (3d = 1 there is always a finite 
fraction of eigenvalues remaining real, as for /i = the eigenvalues of Ip lie along the imaginary y-axis. 
For (3d = 4 a gap opens up on the y-axis, depleting the eigenvalues. If we were to look along the real 
x-axis on the same scale - note the difference of scale in x and y direction of the picture - we would see 
the same depletion here (see e.g. fig. [5] in section d]). Our aim is to compute all complex eigenvalue 
correlation functions analytically in the vicinity of the origin and compare them to Lattice simulations 
where possible. We have succeeded so far for the 2 classes (3d = 2 and 4. 

The addition of the /U-part in the MM in tabled] is not unique. Instead of adding /i times unity Ij, 
as was first introduced by Stephanov [12], the symmetries also allow to add a second random matrix 
in the /i-part, Ij — > of the same symmetry as in the corresponding class (3d |50j : 




This choice of having 2 random matrices (two-MM) corresponds to a non-diagonal chemical potential 
term in random matrix space. Introduced by Osborn [50] for (3d = 2 and the author [51] for (3d = 4 
this apparent complication has the virtue to allow for a Schur decomposition of -0(/i), leading to 
complex eigenvalue model. In the former model with = Ij [32] and only 1 random matrix (one-MM) 
eq. (|2.ip this is not possible. The two models agree for the partition function and density as we will 
see later, illustrating the concept of MM universality. The same universality is expected to hold also 
for a non-Gaussian distribution in eq. (|2.2p but is much harder to prove, for results on non-chiral 
MMs see [52]. 

With all ingredients given so far we can write down our MM partition function with Nf quark 
flavours of masses rrif and chemical potentials /ij in a fixed topological sector v as follows 



The fermions integrated out in the field theory lead to the determinants of the Dirac operator. Here 
we have chosen \I> = 1] , and fixed the inverse variance to be ./V in order to have a compact support in 
the large- iV limit. 

These chiral complex ensembles introduced in [42] for (3d = 2 and [47] for (3d = 1,4 are to 
be contrasted with the 3 Ginibre ensembles for complex eigenvalues introduced in [53] . In [53] no 
determinants are present and the eigenvalues of the matrix $ are considered, instead of those of the 
chiral matrix Ip{p) in eqs. (12. lh or (12. 4p . 

2.2 Relation to field theory 

Apart from using global symmetries in our MM approach it appears so far that we have made the 
approximation of static gauge fields with Gaussian fluctuations, in order to arrive at our MM eq. (|2.5p . 
In reality our approximation is that of static meson fields, and since the first introduction of MMs 
in 1993 [1] this approximation has been understood much more precisely. In the large-iV limit MMs 
become equivalent to a theory of mesons in the limit where their zero-momentum modes completely 
dominate the partition function. This is the e-regime (or microscopic domain) of as explained by 
Gasser and Leutwyler [18], and we will point out this close relation now. Unfortunately, very little is 





* y 6 M/C/real M for (3 D = 1/2/4 



(2.4) 




(2.5) 
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known about the corresponding group integrals over the Goldstone manifold in the classes (3d = 1 and 
4: even at \i = analytic results only exist for degenerate masses for more than Nf = 2 flavours |17j . 
Therefore we will restrict our discussion in this subsection to the (3d = 2 symmetry class including 
QCD, where most explicit results are known. 

The starting point is our Gaussian MM partition function eq. (|2.5p . In a first step we rewrite the 
determinants in terms of Grassmann vectors x, using detLA] = J d\dx^ ex P[ — TrxAx ]• The random 
matrix $ can now be integrated out, being quadratic in the exponent. This leads to a fermionic 
integral with quartic terms in the exponents. The introduction of an auxiliary complex matrix Q of 
size Nf x Nj, also called Hubbard-Stratonovic transformation, makes the fermions x quadratic again, 
and we can integrate them out to arrive at 

ZW") ~ f dQdQUet&Ydet[QiQ-QiBQ^ l B] N e- N ^® f - M M- M \ (2.6) 

after shifting Q — > Q — M for real masses mt. Details can be found e.g. in [54J. Here we have 
abbreviated by M = diag(mi, . . . , mjv"/) the diagonal mass matrix and by B =diag(/ii, . . . , un s ) the 
Baryon charge matrix, the diagonal matrix of chemical potentials. This representation also called 
sigma model is exact, so far no approximation has been made. Notice the duality with eq. (12. 5h . 
where N and Nf have changed role: we now have a determinant to the power N susceptible to a 
saddle point approximation. In particular after introducing also a temperature T the discontinuities 
of the partition function leading to its phase diagram can be analysed in this way when making further 
assumptions about the saddle point, as will be discussed in the next section. 

The partition function can be computed without further assumptions, leading to the chiral La- 
grangian in the e-regime. If we choosqj to scale the chemical potentials keeping 2N/j/j to be fixed in 
the large- iV (volume) limit we obtain 

Z^ N f^ ~ e -ATrA/ 2 



x exp 



J dQdet[Q^e- NTr ^Q 

NTi (M(Q + Qt) - Q~ x BQF~ x B\ + iVTrln(QtQ) + OQ./N)] , (2-7) 

where we have expanded the logarithm. Next we parametrise Q = UR with U £ U(Nf) unitary and 
R Hermitian positive definite matrices. In the saddle point limit we obtain R ~ 1]^ (after going to 
eigenvalues, see [54J) so that the Gaussian term in QQ^ and the logarithm term disappear. Thus we 
can match with the following group integral of e^PT in the fixed sector of topological charge is, up to 
constant terms ~ expfconst • fj?], 



zI"pt= [ dU det[UY ^\\{qq)VT,M{U + U- 1 ) - ^^ 2 Tr[C7, B) [U~ x , B] 



(2i 



It includes the pion decay constant F n , and the chiral condensate (qq). The second, constant term in 
the commutator, Tr[U, B][U~ X , B] = 2Tr{UBlJ- 1 B - B 2 ), can be included in eq. (12. 7ft by multiplying 
with an appropriate constant. This leads us to the following identification of MM and physical, 
dimensionful parameters: 

2Nm f — ► V(qq)m f (2.9) 

In principle we could have introduced two parameters in the MM through its variance, e~ N ^ Tr*^ ; 
and by rescaling fj,f — > const- [if . At fi = the MM predictions measured in units of {qq) are parameter 

A different choice is possible, keeping fXf fixed at iV — > oo. This leads to a different result, see [54] for details. 
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free. At fj, ^ the ratio of the two parameters has to be determined, e.g. by a fit to lattice data. This 
is is a virtue and not a drawback, as such a fit determines F n without having to look at finite-volume 
corrections (55[ [56] . 

Let us explain how the group integral eq. (|2.8p is obtained in the e-regime of xPT, following 
[18, 16J for /i = and [20j for fi ^ 0. To leading order the chiral Lagrangian is given by 

£xPT = ^fV a U{x)V a U{x)^ - ±(qq)(Me^U(x) + (Me™tU(x)?) (2.10) 

Here we have included a topological 0-term in the QCD gauge action S gauge : +i@ J FF = i@v, related 
to the exact Dirac operator zero modes due to the index theorem. The latter produces a factor m'j 
from the Dirac determinant for each flavour, which is why mass and G-term always group as shown. 
The space-time dependent field lives in the coset space U(x) E SU(Nf). The covariant derivative 
V a U = d a ll — n(B a U + UB a ) introduces a coupling to the chemical potentials via the Baryon charge 
matrix B a = (B, 0). In the chiral limit mj —* the zero momentum modes lead to a divergence in the 
propagator. It can be cured by splitting off these zero-modes U and treating them non-perturbatively: 
U(x) = Ue^ x \ where U 6 SU{N f ) is now a constant matrix. The propagating modes <p(x) have 
to be expanded in the algebra, that is in terms of the Pauli matrices for Nf = 2 for example. In 
the e-regime of xPT in a finite volume box V = L 4 [18j the pion mass is counted as m n ~ 1/L 2 (in 
contrast to the p-regime where ~ 1/L): 

m n ,fJL <C 1/L -C A . (2.11) 

The right inequality implies the validity of the chiral Lagrangian, with A being the scale of the lowest 
non-Goldstone modes. The left inequalities imply that the path integral over U(x) factorises^l into a 
constant SU(Nf) group integral times a Gaussian integral over the fluctuations of (p(x). The latter 
gives only a constant that is omitted. In a last step eq. (|2.8|) is obtained by fixing the sector of 
topological charge in the sum Zqcd = Yl°^oo e i&v Z( N f' u ': the inverse Fourier transformation leads to 
the prefactor detfC/] 1 ' and promotes the integral from SU(Nf) to U(Nf). Let us emphasise that in the 
e-limit the leading order Lagrangian eq. (12.101) becomes exact. 

What we also learn at this point is when the MM approximation breaks down: it happens when the 
zero-momentum Goldstone modes cease to dominate and propagating modes start to contribute. In 
analogy to condensed matter this scale was called Thouless energy, and it is given by E c / (SV) ~ F^W 
in dimensionless units at /U = [22]. Clearly \i ^ will influence the mixing and it was found in |57j 
that the Thouless energy increases with \x on the lattice. 

After pointing out the equivalence between MMs and the e^PT on the level of partition functions 
- a result which is new for more than 2 chemical potentials and was independently derived in [58J - we 
give some explicit results in which the partition function can be further calculated in terms of Bessel 
functions. For m + n = Nf quarks of opposite isospin chemical potential B = fx diag(l| m , — l] n ) with 
rescaled masses {rjj± = V {qq)nij±} respectively, we obtain for the group integral eq. (|2.8p . as well as 
for the MMs eqs. (JZI]) and ([23]) [59] [53]: 



7 (N f ,u) 1 , 

^ PT ~ A m ({ v l })A n ({e }) l<ff<n 



(2.12) 



Without loss of generality we have chosen n > m here. With A m ({r/ 2 }) = Hfc>z( r ?fc — r n) we abbreviate 
the Vandermonde determinant of squared arguments. Here we have introduced the notation 



f 1 dt t e- 2VF ^ 2t2 I u {tr, + )I u (tr,_) , (2.13) 
J o 



6 The terms mixing propagating and zero modes are boundary contributions [56] . 
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constituting the first m rows. The remaining rows contain lj) k \r]) with increasing k, denoting the 
kth derivative of the modified I-Bessel function with respect to its argument. Note that the partition 
function eq. f)2. 12j) is real only if all masses are real (or purely imaginary). 

When all quarks have the same chemical potential fj,f = fj,, we recover the familiar result [16] 

Z)L T ' ~ det 

^ PT A Nf {{rj 2 _}) l<f,g<N f L 



a / I 1 /(/- 1) (%-)l , (2-14) 



which is totally independent of /x. This comes as no surprise since pions don't carry Baryon charge: 
[U, B ~ 1] ] = 0. In particular this implies that the Leutwyler-Smilga sum rules [16] continue to hold, 
despite the eigenvalues being now complex. The simplest example is given by 

(2.15) 




A{v + N f ) ' 

for Nf massless flavours, for massive sum rules we refer to |60j . Here X' denotes the sum over ^ 
only. For quarks including isospin partners sum rules depend on p, as has been worked out explicitly 
for imaginary /i [61]. Note that the group integral eq. (|2.8p equally holds for imaginary chemical 
potential by the simple rotation /i — > i/i, flipping the sign of /i 2 . On the level of the initial MM such 
a rotation is highly nontrivial, and for the corresponding MM and its solution we refer to section [6l 

Let us make an important remark. The computation of the group integral eq. (12, 8D is relevant 
even for a theory of only pions: despite the fact that its partition function is /x- independent, the 
generating functional for its spectral density P^(z) in the replica approach [121 [59] does depend on /x 
through a pair of conjugate quarks of masses z and z* . This statement directly relates to the failure 
of the quenched approximation as explained in [32], and will become more transparent in subsection 
O later. 

In addition to the above also partition functions with one or more bosonic quarks have to be 
computed in the Toda lattice approach [59J, giving rise to inverse powers of determinants. In the 
simplest case of 1 pair of conjugated bosonic quarks, denoted by negative Nf = —2, the partition 
functions is reading [59] 



A e X PT - J det[Q] 2 e e K y J 

(2.16) 

for zero topology v = 0. The integration goes over 2x2 positive definite Hermitian matrices and has 
to be regularised. Only the leading singular term is given here, which can be reproduced from MMs 
as well [63]. The analytic behaviour and //-dependence of such bosonic integrals is strikingly different 
from fermionic partition functions, and we refer to [M] for details. 

General partition functions with both fermions and bosons can be also used directly as generating 
functionals for resolvents in the supersymmetric approach. This has been used at /i = in |19|, I20|. in 
order to derive the spectral density of the Dirac operator directly from e^PT, without using MMs. In 
the following sections we will use a different approach that does not rely on bosonic or mixed partition 
functions. 
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3 Phenomenological Application: Phase Diagrams of QCD 



In this section we will derive a qualitative phase diagram for QCD using the MM partition function in 
terms of the sigma model from the last section. These results are called phenomenological as the link 
to the chiral Lagrangian is lost close to the phase transition, where an d in particular the epsilon 
regime breaks down. Let us point out however, that for the two other classes, SU(2) and adjoint 
QCD, such an effective description can be used along the /i-axis to describe a phase transition due to 
diquark condensation. We refer e.g. to [31] for further details. 

Coming back to QCD let us write down the effective Lagrangian or potential, following closely 
|65j. The effect of the lowest Matsubara frequency ±7rT or Temperature is introduced as in [65J, 

Th( rp, ( $ + (fl + iTa 3 )^ \ 

^{-& + (» + iTa 3 n J' (3 - 1} 

Because both the Euclidean time derivative and the //-term are proportional to 70, (d T + (J,)jq, they 
are replaced by a constant matrix times an off-diagonal unity matrix. In addition to the first proposal 
J49j here half of the Ip eigenvalues have positive and half have negative frequencies, through the 
<73-term. Despite the dimensional reduction through the discretisation of the time direction the the 
global symmetries of Ij) are taken to be the same as in 4 dimensions. Note that the above T-term 
alone would preserve the anti-Hermiticity. There are ways to include also higher frequencies by T — > 
diag(Ti, . . . ,T/v) as suggested in [66], but we will pursue the simplest approximation [65]. Also we 
only keep one single, real fi here for all quarks. All steps go through as in the previous section until 
eq. (12. 6p . In terms of the auxiliary field Q we obtain the effective Lagrangian 

£(Q,{m},/i,T) = N{qq)QQ^-^ln[(Q+M)(Ql+M)-(fi+iT)][(Q+M)(Q^+M)-(Li-iT)]. (3.2) 

The topological charge v of order 1 can be neglected in the large- N limit. Making the further assump- 
tion that the saddle point solution is proportional to unity, Q = 4> lj , leads to the following situation. 
Despite C not being polynomial the resulting saddle point equation S^C = is a polynomial equation 
of 5th order, just like following from a Landau-Ginsburg ansatz with a sextic potential. At zero mass 
one obtains [65] 

<t> 4 - 2 (V - T 2 + <f) 2 + (/i 2 + T 2 ) 2 + /j 2 —T 2 = , (3.3) 

leading to a second order line starting at /i = and T c to meet a first order line starting from T = 
and \i\ in a tri-critical point (T3, ^3) = (5V \/2 + 1, \\J y/2 — 1) [65]. For equal nonzero mass M = ml] 
the second order line becomes a cross over as shown in the full 3-dimensional phase diagram fig. [2j 

Using the input of T c = 160 MeV and \i\ = 1200 MeV physical scales can be reintroduced, leading 
to a prediction for the location of a tricritical point. 

Some remarks are in order. The same phase diagram has been obtained from several other models 
as reviewed in [68], as well as partly from Lattice simulations along the T-axis and close to T c following 
the transition (crossover) line (see \29\ [32] for references). However, there is no flavour dependence 
in this MM prediction as the dependence of the sigma model on Nf is too weak: the phase diagram 
looks the same also for any number Nf > 2, given there is no other saddle point than Q = On 
the other hand for Nf = 3 flavours we know from Lattice studies that for zero [i the transition at 
T c depends on the masses, e.g. for equal mass it becomes first order. We refer to [69J for a recent 
review. As a further restriction (qq) is a fixed, Nf independent input parameter in our MM. Because 
of the difficulties of performing Lattice simulations with a complex action, the presence or absence of 
a tri-critical point is still a subject of debate today (see [32] for a most recent critical review). 
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Figure 2: Phase diagram of QCD for Nf = 2 light flavours, from [65J. 



The MM can be refined by including different chemical potentials for the quarks, that is baryon 
(hb) and isospin (///) chemical potential. This analysis has been carried out in [67] for Nf = 2 
including a very detailed analysis of the phase diagram. We only cite the most prominent feature, 
that is a doubling of the phase transition lines for fixed fij and mass m as shown in fig. [3l This can 




Figure 3: Phase diagram of QCD for Ay = 2 light flavours in the \ib plane for fixed \xi and mass m, 
from [67J. All quantities are in units of the inverse variance G. 

be easily understood as for fj,\ ^ fi2 the chiral restoration for the up and down quark condensates 
(uu) and (dd) is different. A similar effect can be expected if the masses are split, m u / m^, at equal 
chemical potential \x\ = The same split of one first order transition into two was found analytical 
from MMs at T = [54j . as can be derived from eq. (|2.7p of the previous section. It leads to the 
fact that the iVy-flavour partition function can be written as a determinant of single flavour partition 
functions. In consequence the Ay-flavour partition function undergoes a transition whenever one of 
its building blocks, a single flavour partition function, has a phase transition (singular derivative). 
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4 Exact Limit: QCD Dirac Operator Spectrum 



In this section we consider our MM as an exact limit of QCD, in the sense that it coincides with the 
e-regime of as has been pointed out in subsection 12.21 This assumes that xPT as an effective 
theory itself can be consistently derived from QCD. In all the following we will consider temperature 
T = only. 



4.1 Chiral symmetry and the Dirac spectrum 

To start we recall the relation between the spectral density of the Dirac operator and spontaneous 
chiral symmetry breaking from field theory, both for zero and non-zero fi. Throughout this subsection 
no reference to MM results is made, apart from figs. H] and 

Let us begin by defining the spectral density of Dirac operator eigenvalues eq. (|2.3p by 



\ * / QCD 

where the average is taken over the gauge fields including the Nf Dirac determinants. In our conven- 
tions the eigenvalues z^ are real at [i = 0. In the same way higher order correlation functions can 
be defined, such as density-density correlations, etc. . In a finite volume on the Lattice there is only 
a finite number of discrete Dirac eigenvalues. When computing integrals at limy^oo with such an 
operator, such as sum-rules eq. (j2.15fl . special care has to be taken for UV divergencies, and we refer 
to [70] for this issue at \i = 0. Putting QCD on the Lattice always provides a natural regularisation 
here. Because of chiral symmetry 

{(0 + /i 7o ),75} = , (4.2) 

eigenvalues of (lp + fi^o) come in pairs Hz^ as for each eigenfunction tp^ of a non-zero eigenvalue there 
is also an eigenfunction j^ipk with eigenvalues —iz^, except for zi~ = 0. The number of zero eigenvalues 
relates to the topological charge, v = \til — tir\. 
Because of this we can write for the resolvent 

(4.3) 

where the second sum E' goes only over the non-zero eigenvalues with —tt/2 <Arg(z&) < vr/2, the 
half plane C+. Here m is an auxiliary mass. The partition function of course depends also on the Nf 
masses mf, see also the discussion after eq. (|4.33p . The chiral condensate is then defined as 

(qq) = — lim lim S(m) , (4-4) 

where the order of limits is important. Using the electrostatic analogy £(m) is the electric field for 
charges izt on the imaginary axis for [i = 0, or on a stripe parallel to it for /i ^ 0. The fact that chiral 
symmetry is broken trough the chiral condensate (qq) corresponds to a jump of £(m) along this axis, 
see fig. H] below. 

The relation between resolvent and spectral density however is different for Zk real or complex, as 
summarised in table [3l For zero chemical potential we have 

y limE(m)=/ d^-j—, (4.5) 
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n = o 




S(m) 


lim n _ ±d m (det(Ip + m) n ) QC D 


lim n ^ ^ m (det + m\ 2n ) QCD 


limy^oo E(m) 






Pjf z = x + iy) 


lim e ^o gSiP(*m + e) - E(im - e)]\ m=x 


dm*^(m)\ m=z 



Table 3: The resolvent and spectral density for zero and non-zero fi. Here m is an auxiliary mass, the 
iVj-dependence of pm and S(m) has been suppressed. 



with x 6 K.+, assuming that (u 2 ) ~ V is suppressed. This leads to the Banks-Casher [7T] relation in 
the chiral limit 

l<OT>l = £/ty(0) , (4-6) 

where we have used the representation lim m ^o Ji+^i = ^^(x) of the 1-dimensional delta-function. 
On the right hand side of eq. (|4.6p we have the macroscopic density of Dirac eigenvalues at zero 
momentum. Its slope has also been computed in [48J, p'^{x = 0) = 1 j^] F 4 — jv)fo) + ^ ? depending 
on the symmetry class through the Dyson index. Because of eq. (|4.6p the smallest eigenvalues are 
spaced like 1/V (in contrast to free eigenvalues ~ l/V*) and we define a microscopic density, both 
for \x = and \x ^ 0, 

PS(£) = M£ = z/(qq)V)/(qq)V . (4.7) 



This microscopic density can be obtained from eq. (|4.5p by inverting the integral equation - if p(x) is 
independent of the valence or source quark mass m- as the discontinuity of the resolvent: 

Pj j(x) = lim — \E(ix + e) - + e)] . (4.8) 



On the other hand for non-zero ji we have 



+ 



2m 



lim E(m) = / d 2 z pj(z)-5 = . (4.9) 



Here the chiral limit m — > does noi produce the 2-dimensional delta-function tt5 2 (z) = d z *~. The 
density - resolvent relation in the complex plane instead reads: 

Pjfz) = \d z *Y,{z) (4.10) 

that is outside the support the resolvent is holomorphic and inside it depends on the complex conjugate 
z*. For more details we refer to [72J. When adding a delta- function <5 1 (y) to eq. (|4.8p we could also 
write it in the form of eq. (|4.10p . 

So far we have only defined the microscopic density and the resolvent, and related the two. The 
question remains of how to compute them from field theory. We give two possibilities, replicas and 
super symmetry, where again we have to distinguish // = and fi 7^ 0. Here extra fermionic and 
bosonic valance or source quarks are added to generate the resolvent. 
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Figure 4: The microscopic spectral density p s f ' v (£) at large volume, cut perpendicular to the real 
axis and rotated by tt/2: the real part of the complex valued, unquenched density from MMs for 
Nf = 1 flavour with mass m(qq)V = 80 (upper plots) vs. the quenched density from e^PT (lower 
plots), both at rescaled fiF^y/V = 8. The resolvent T,(m(qq)V) obtained numerically is shown for 
comparison. All plots are from the second ref. of [73], with S = (qq). 



In the supersymmetric approach at p = we can write 



E(m) = d, 



mi 



/det(# + mi) 
\det(-0 + m 2 ) 



QCD 



(4.11) 



m\=mi 



using the fact that det(0 + m) = exp[Trlog(0 + m)]. Here mi and m,2 are 1 extra fermionic and 
1 bosonic quark, respectively. As long as m,2 has a small imaginary part, the generating functional 
containing the ratio of determinants in addition to the Nf flavours is well defined. In the limit of 
exPT we obtain an integral over the super Riemannian manifold Gl(Nf + 1|1) generalising eq. (|2.8p . 
and we refer to [19] for details. 

For [i the inverse determinant can no longer be regularised by an imaginary part, as the 
spectrum of lft(n) now extends into the complex plane. Here the so-called Hermitisation technique is 
used as follows. Defining 

'det[(0(/i) + mi)(0Gu) + m x )t + k 2 } \ 

I QCD 



det[(#(/i) + m 2 )(Ip(fi) + m 2 )t + k 2 } 



(4.12) 



we can obtain the regularised resolvent by differentiation and setting arguments equal as before, and 
then the density through eq. (j4. 10H : 



P 



1 



lim cL 



7T K^O 



lim d m2 Z( N f +2 ^\m 1 ,m 2 ; K ) 



ni2— >m\ 



(4.13) 



mi- 
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For a detailed calculation in a non-chiral MM framework we refer to [75J. 

When using replicas the logarithm of the determinant is generated by inserting n fermionic deter- 
minants of degenerate mass (replicas) into the partition function, leading to the resolvent as follows: 

S(m) = lim -d m (det(# + m) n ) nrn . (4.14) 

We will not comment here on the subtlety of analytically continuing n — ► 0, as we will not use replicas 
in our computations. Eq. (|4.14p is only valid for real eigenvalues at [i = 0. For \i ^ the appropriate 
replica limit includes complex conjugated quarks 

E(m) = lim -d m (det[(lp(u) + m) n (IpU) + m) tn A . (4.15) 

n— >0 n \ I QCD 

This was pointed out by Stephanov in [32] to explain the failure of the quenched approximation in LGT 
at \x 0. Using MM eq. (|2.5p he showed that the critical potential is proportional to the pion mass, 
\x c ~ m w , due to the presence of conjugated quarks in the correct replica limit eq. (|4.15p . This has to 
be compared with a third of the lightest Baryon mass, p? c ~ m^/3, expected for unquenched QCD. 
Replicas for non-hermitian MMs and e%PT have been introduced in |73[ [59] to compute correlation 
functions, using a relation to integrable hierarchies to make the replica limit exact. 



X(m(qq)V) 




10 -10 2 4 6 8 10 12 14 



Figure 5: The same cut as in fig. 2] for the quenched density in the (3d = 4 or adjoint QCD class from 
MMs at rescaled fiF n y/V = 2.5 and topology v = 0. The corresponding resolvent is shown for positive 
arguments only. 

We close this section by giving cuts through the microscopic densities in the large volume limit as 
obtained from MMs in the following sections. This is to illustrate the different symmetry classes and 
the effect of unquenching. Parallel to that the resulting resolvent is shown. 

As was pointed out in [73] the strong oscillations in the complex valued unquenched spectral density 
are made responsible for chiral symmetry breaking, with a discontinuous resolvent. For comparison the 
quenched density and its continuous resolvent are shown as well in fig. [H Here the chiral condensate 
is vanishing in the chiral limit and thus not responsible for symmetry breaking. It is the diquark 
condensate breaking chiral symmetry in this case. 

We also give a cut through the quenched density for adjoint QCE0. The resulting resolvent is again 
smooth as in quenched QCD above. This implies that the chiral condensate (qq) is rotated to zero, 
with chiral symmetry being broken through a diquark condensate. 

7 Because of slow convergence for the double integral in density eq. (|5 . 6|) we give smaller volumes compared to fig. [4] 
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4.2 Complex Dirac eigenvalue correlation functions from MMs 

In this section we will show how to derive complex eigenvalue correlation functions from MMs using 
orthogonal polynomials in the complex plane. The reason we choose this technique is that it is simple 
and carries over with little technical modifications to symplectic MMs or adjoint QCD in the next 
section. The alternative techniques are replicas using the Toda lattice, which has been successful only 
in the (3d = 2 class, and supersymmetry. They will be mentioned only very briefly. 

Our starting point is the two-MM version of the partition function eq. (|2.5j) based on eq. (12. 4h . 
that was first introduced in [50j (see also [63] for details) 

Nf 

= / n « ( viv * +;* ) e -^<***+«'> (4.i 6 ) 

In order to obtain a complex eigenvalue model a Schur decomposition for the following complex 
matrices is performed, C = $ + /i* = U(X + R)V and D = — $t + fi^f = V^(Y + S)!/^ . Here X and 
Y are diagonal matrices of complex eigenvalues x k and y k . R and S are complex upper triangular 
matrices, U and V are unitary. The Jacobian is proportional to An({z 2 }) = Ylj>k( z j ~ z k)> the 
Vandermonde determinant. Because of the Gaussian measure the matrices R and S decouple and can 
be integrated out. A change of variables to the complex Dirac eigenvalues z k = x k y k of the matrix 
CD is made next, and finally one of the 2 complex eigenvalue sets, {x k } is integrated out. This last 
step is transforming part of the Gaussian measure into a -fT-Bessel function: 

z(Nf,u) =J c fi d2z i n mv M +™/) i%-r +2 ^ ( N{ \y 2) ^ 2 ) e^^ 2+2;2 W({^})i 2 • 

11 (417) 

Here we have omitted all normalisation constants as they will drop out later in correlation functions. 
An approximation to this complex eigenvalue model was first introduced in [76] replacing K v {x) — > 
a/7t/(2x) e~ x by its asymptotic value. The two models coincide at v = ±5. 

The concept of orthogonal polynomials (both on M. and C) to compute eigenvalue correlation 
functions relies on the following properties. Given a weight function w(z,z*) on C such that all 
moments are finite we can define a set of monic polynomials Pk(z) = z k + 0{z k ~ l ): 

d 2 z w(z, z^PkizWz)* = S kl h k (4.18) 

c 

where hu are the squared norms. For real positive weights the coefficients of the polynomials can 
easily be seen to be real, Pi(z)* = P[(z*). Using invariance properties of determinants we can write 
the following identity: 

N-l 

A N ({z 2 })A N ({z* 2 }) = del : [P k ^{ Zl )] det [JV-i(*f)] = \\ h f , Jet [K N (z k ,zt)} (4.19) 

Kk.KN Kk,KN J \<kl<N 

- - - ~ j= - - 

where we have defined the kernel of orthogonal polynomials 

N-l 

K N {z k ,zt) = Y, h J lp ^k)PM) ■ (4.20) 

3=0 

It enjoys the following self contraction and normalisation property: 

d 2 z w(z,z*)K N (z k ,z*)K N (z,z*) = K N (z k ,z*) 

... 

d 2 z k w(z k ,z* k )K N (z k ,zl) = j ^ !j k k = J . (4.21) 
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For such kernels the following Theorem by Dyson (Theorem 5.2.1 in [44]) holds 

[ d 2 z n w(z n ,z* n ) det [K N ( Zi ,z*)) = (N-n + 1) det [K N (z h z*)] . (4.22) 
Choosing the polynomials orthogonal with respect to the weight 

w^\z,z*) = Uy + m}) \z\*>+*K v (^^M 2 ) exp (^A* 2 + ^)) , (4-23) 

in eq. (|4.17p we can apply the Dyson Theorem to compute all complex eigenvalues /c-point correlation 
functions defined as 

/vi i r N N 

R^\ Zl ,...,z k ) ee {N : k)lz[Nf , v) j c II d 2 Zj l[w^\z h z;)\A N ({z^ (4.24) 

k 

= l[w( N f>»\z h zt) i det , [K N ( Zi ,z*)]. (4.25) 
i=i - hJ - 

It expresses the original (N — fc)-fold integral as a determinant of size k x k only. Thus in the large- N 
limit we only have to compute the asymptotic of the kernel Kn(zi, z*). Note that for brevity we have 
suppressed the dependence on the complex conjugated arguments . . . , z% and on the masses in the 
functions R( N f> v \ 

Let us give an example. For the quenched weig ht, w^(z,z*) with Nf = in eq. (|4.23p . the 
orthogonal polynomials are given by Laguerre polynomials in the complex plane [50J 

P^) = H^Cl-M 2 )^^^) • (4-26) 

For a proof of their orthogonality we refer to appendix A in |51j . The corresponding kernel is given 
by 

KMH.4) = -(^iHvf ^ I JJ^L", ( TT7 ) V ( TT7 ) ■ (4.27) 

The orthogonal polynomials for weights including Nf > flavours and their kernel can be expressed 
in terms of the quenched quantities eqs. (|4.26p and (|4.27p using the technique of bi-orthogonal poly- 
nomials [771 150] . or ordinary orthogonal polynomials [78] in the case of a positive definite weight 
(corresponding to Nf > phase quenched flavours). From this all correlation functions follow for 
Nf 0, and for most explicit expressions we refer to |63pl . see also eq. (|4.38p below. The partition 
functions with Nf flavours themselves can also be expressed in terms of the quenched polynomials 
and kernel [75], resulting into eq. (|2.12p in the large- N limit. 

As pointed out in subsection 12.21 the MM only becomes equal to the chiral Lagrangian in the 
following large- N limit. We rescale complex eigenvalues, masses and chemical potential as in eq. (|2.9p 

= 2Nzj , r] f = 2Nm f , and a 2 = 2Nfi 2 , (4.28) 

which is denoted by the microscopic limit at weak non-Hermiticity. The concept of weak non- 
Hermiticity introduced in [79] for non-chiral MMs implies that /i 2 vanishes at a rate 1/N. There 

8 Note the different convention there, rotating eigenvalues by 7r/2: iz —* z. 
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exists a different iimit cailed strong non-Hermiticity with a different scaling for the eigenvalues, keep- 
ing p? = 0(1). This leads out of the domain of the chiral Lagrangian, and we refer to [5l] for 
details. 

The microscopic spectral density in the large- iV limit eq. ([4.28P is defined as follows 

= lim -LRWfrifz^A-) . (4.29) 

All higher order correlations are rescaled accordingly. The quenched microscopic density obtained 
from eq. (|4.27|) thus reads 

^ M(e) = (S) ^ J! dt 1 J M)JM*) ■ (4-30) 

In the Hermitian limit a — > it reduces to the density of real eigenvalues p^' u \£, x ) P] times 5 1 (9m(^)) 

P^HCx) ~ |^| J dt t J v {ti x ) 2 = ^{J v {jU? ~ J v -i&)J v +\(S*)) , (4-31) 

where 3ie(£) = £ x denotes the real part. Fig. [6]compares the two functions for different values of v and 
a. For small a the density in the complex plane eq. (|4.3U|) is very similar to the real one eq. (|4,31|) . 
times a Gaussian decay in imaginary direction. Introducing exact zero eigenvalues by increasing v 
pushes the density away from the origin through level repulsion. For increasing a the density quickly 
becomes a constant along a strip as from mean field, washing out the oscillations. 

For the microscopic partition function with Nj flavours we obtain eq. (|2.14p which is /i-independent 
for B = 1| as it should. We can now check that the same sum rule eq. (|2.15p indeed follows from the 
quenched density eq. (|4.30|) : 

4^ Jc C 

which we have confirmed numerically up to v = 5 for two different values of a < 1. Thus despite that 
for n ^ the eigenvalues spread into the complex plane the sum rule remains unchanged. An analytic 
check, in particular of the unquenched, massive sum rules [60] using the complex valued densities for 
Nf > 1 given below would be highly desirable as well. For fx = the Leutwyler-Smilga sum rule eq. 
(|2.15|) was shown analytically to follow from the real density eq. (|4.31|) |80j. 

Before turning to more flavours let us comment on other methods of deriving complex eigenvalue 
correlations and their relation to field theory. As was pointed out in subsect. 14.11 the spectral density 
can also be derived using replicas [59] 

pf f '"\0 = lim— d^\og[Z^f +2n ' v \{T,},{^*})] , (4.33) 

where the first derivative generates the resolvent. Using replicas for complex eigenvalues (see table [3]) 
we have introduced n paird3 of complex conjugate quarks of masses z£ and i^* respectively, having 
opposite (j,. Thus the corresponding partition functions eq. (|2.12p now depend on li. It was shown in 
|59] that such partition functions obey the Toda Lattice equation 

(4.34) 



9 I am indebted to Leonid Shifrin for his help with Mathematica. 

10 Usually a different notation is used, counting only the number of complex conjugated pairs n in the superscript. 
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Figure 6: The quenched microscopic density eq. (|4.30p at a = 0.4 for v = (top left) and v = 2 (top 
right), and at a = 3.28 for v = (bottom right). For comparison we also show the density of real 
eigenvalues eq. (14.31D (bottom left) for v = (left curve) and v = 2 (right curve), respectively. Note 
that in the complex plane the normalisation becomes a delta-function in the Hermitian limit a —* 0. 



leading to 

^ ^ = 5 K W^HM)) 2 ' ( } 

Here the minus 2 means that this partition function contains 2 conjugate bosons instead of fermions. 
In [63] the corresponding bosonic group integral given by eq. (|2.16|) was computed for Nf = 0, thus 
deriving the quenched density from eq. (|4.35j) entirely from e^PT, without using knowledge from 
MMs. For Nj > 1 this result is still lacking, and using MMs we know that the corresponding group 
integral will factorise: Z^f' 2 ^ ( {rj}, {£, £*}) ~ ^ N f' v \{q})Z^ 2 ' v \{i,i*}). The bosonic partition 
function has to be regularised and is proportional to the weight function eq. (|4.23p to leading order. 
Using this as an input the equivalence of the Toda lattice and the orthogonal polynomials approach 
in MMs has been shown in [63] , see eq. f|4.38[) . 

Another method to compute the spectral density is supersymmetry. Here the additional inverse 
determinants are introduced as Gaussian integrals over bosonic Grassmann variables, see the discussion 
after eq. (|4.11|) . and then treated along the lines similar to subsect. 12.21 We refer to [81] where the 
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equivalence of the supersymmetric and the Toda lattice approach was shown. 

After comparing different techniques of solving MMs we would like to comment on the concept 
of universality. Sometimes the fact that instead of using a field theory such as e^PT to compute 
correlation functions one can use a simpler, chiral random MM, is called universality. As we have seen 
in sect. 12.21 these two are equivalent on the level of partition functions. 

Here, we would like speak about MM universality, meaning that in the large- A limit different 
MMs, Gaussian or non-Gaussian, containing 1 or 2 random matrices, lead to the very same eigenvalue 
correlation functions. What is know about universality for MM with complex eigenvalues? So far 
we know that the Gaussian one-MM eq. (|2.5p and the Gaussian two-MM eq. (|4.16p have the same 
partition functions for any number of Nf fermions with chemical potentials ±/i, as well as the same 
quenched spectral density. Higher order and unquenched correlation functions have been obtained 
only in the second model. 

At present little is known about the large-A limit of non-Gaussian chiral MMs with complex 
eigenvalues. The determinantal |76[ l50l [63] or Pfafhan structure [51j. [82] of the correlations functions 
remain the same for non-Gaussian weights at finite- A, and from heuristic universality results [52] for 
the non-chiral complex ensemble at (3d = 2 we expect that they are also universal. 

We finish this subsection by giving the results for more flavours. With the definition 

^(6,6)= / dtte- 2aH2 JMi)JM2) , (4.36) 
J o 

the unquenched density for Nf = 1 reads [631 EQ] 

(.37) 

where we have factored out the quenched density. Note that the support of the quenched and un- 
quenched density is the same. It can easily be seen that the density p^ f is no longer real and 
positive. For general Nf the structure remains the same : the factor multiplying the quenched density 
becomes a determinant of size Nf + 1 [63\ 150] : 



|e| K » x& ) e ^ n^' + f ) — Z (v) (M) — • ( 438 ) 



where the additional pair "+2" has masses i£ and The weight function originates from the bosonic 
partition function in the Toda lattice picture, and we refer to [63J for more details. It is not the explicit 
mass factor in front that makes the densities complex as it gets cancelled by the Vandermonde from 
eq. (|2.12p . We recall that the partition functions themselves become complex valued for complex mass 
arguments £. 

For illustration we give the spectral density for Nf = 2 in different theories. In the upper plots figd 
real and imaginary part start displaying oscillation for moderate values of rescaled chemical potential. 
The spectral densities for phase quenched can also be computed using orthogonal polynomials, and 
wc refer to [63j for details. For two phase quenched flavours with | det + m)\ 2 the density 
with the very same parameter values is given in fig. [7] by the lower left plot. The oscillations have 
disappeared, with the real density simply vanishing at the locations of the mass. Note the difference 
in scale, the height of the supporting strip is hardly visible in the upper left plot. We also show the 
same unquenched density in the (3d = 4 or adjoint class, from section [5j Apart from the additional 
repulsion from the axes it looks very similar to phase quenched QCD. 
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Figure 7: Comparison of densities for different theories in the first quadrant of C (rotated by 7r/2). 
Upper plots: the complex spectral density for unquenched QCD with Nf = 2, $te(p(z)) (left) and 
( ^sm(p(z)) (right), from |63j . Lower plots: phase quenched QCD (left) and adjoint QCD (right), both 
with Nf = 2. All plots are at the same values a = 2, mass r\ = 5. 

4.3 Comparison to quenched Lattice simulations 

In this section we compare the predictions for complex eigenvalue correlation functions from the 
previous section to Lattice simulation. Because of the sign problem only quenched or phase quenched 
simulations can be compared with so far, and we will focus on the former. Phase quenched simulations 
with two flavours have been performed (see e.g. |83j) but the Dirac spectrum has not yet been 
compared to MMs. 

The first comparison with the spectral density of quenched QCD on the lattice was performed in 
|84j for 3 different volumes 4 4 , 6 4 and 10 4 using staggered fermions at gauge coupling = Q/g 2 = 5.0 
for chemical potentials varying from p = 1CT 3 - 0.2. The reason for choosing a rather strong coupling 
was to make the window in which MMs apply large enough for the small lattices studied. This assumed 
that the concept of a Thouless energy [22], the maximal eigenvalue up to which MMs and exPT agree, 
can be generalised to fi ^ 0. Evidence for this picture has been provided later from Lattice data in 
[57]. 



23 




Figure 8: Comparison of the quenched density eq. (|4.30p at a = 0.2 (full line) and Lattice data 
(histograms) for V = 6 4 and \i = 0.006 (upper plots), from |84j . The left curve is cut along the real- 
axis, and the right curve is on the first maximum of the left curve, cut perpendicular to the real-axis. 
Fitted integrated density (lower plots) eq. (|4.39p y s- small Dirac eigenvalues from Lattice data [57] 
for V = 6 4 and \i = 0.1 (left) and [i = 0.2 (right). The two curves are indistinguishable. 

Instead of comparing to the 3-dimensional plots of fig. [6] we show cuts compared to data. For 
small a shown in fig. [8] we cut along the real axis and parallel to the imaginary axis on the first 
maximum. The agreement between the data and eq. (|4.30p is excellent. In the plot the asymptotic 
form of the K-Bessel function in eq. (|4.30p was used, as predicted from [76]. The exact expression eq. 
(|4.30p was only found later in [59]. However, at the fitted value of a = 0.200(2) the two curves are 
not distinguishable from the data (see fig. 1 in the second of ref. |59j). 

For all other plots we use the exact density, eq. (|4.30p . The lower plots of fig. [H]from [57] show 
reanalysed and extended data from [84] . At larger values of a the density at the origin becomes 
asymptotically rotational invariant, see also fig. [9] Therefore the integrated density eq. (|4.30p 

l(r,0)= sds d<t> ps(t = se**) . (4.39) 
Jo Jo 

is compared to the data, finding again excellent agreement. Because of the rotational invariance only 
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the ratio Ti/F^ can be determined, and we refer to [57] for details. 

Recently progress has been made concerning chiral fermions with chemical potential. In [85] it 
has been shown how to introduce [i ^ into the overlap operator [86] using a generalisation of the 
sign function in the complex plane. This opens up the possibility to compare to different topological 
sectors, and successful comparisons to the density eq. f|4.30j) have been made for v = and v = 1 
shown in fig. Here the same cuts as in fig. [5] have been made (real and imaginary axis have been 
interchanged rotating by vr/2), finding very good agreement. For the simulation and details of the 
complex overlap operator, in particular concerning the projection from the complex Ginsparg- Wilson 
circle to the complex plane we refer to |85j . 



H = 0.1 




Figure 9: Cuts through the quenched density eq. (I4.30|) for topological sectors v = and 1, vs. Lattice 
data from [85] for V = 4 4 and \x = 0.1 (upper plots) and \i = 1.0 (lower plots), (x = 5fte£ and y = 9m( 
are interchanged w.r.t. fig. E|). In the left plots the vertical bar indicates the fit range. For the largest 
value [i = 1.0 the density is almost rotational invariant, as indicated by the lower right plot fitted 
to the asymptotic density ps(ICl) = l£| 2 /2ita 2 K v {\£\ 2 / Aa 2 )I v (\£\ 2 /4a 2 ). In this case only the ratio 
{qq)/F^ can be fitted. 



The results presented here are not the first comparisons between non-Hermitian MMs and Lattice 
data. In [87] quenched QCD lattice data were compared for a large range of \i to predictions of the 
non-chiral (3d = 2 Ginibre ensemble for the spacing between eigenvalues in the complex plane. These 
predictions are only valid away from the origin where chiral symmetry is irrelevant, in the so-called 
bulk region. Here no relation to an effective theory as the chiral Lagrangian is known. A transition 
between the Wigner surmise for fi = 0, the Ginibre, and the Poisson distribution in the complex plane 
was found [87j . Unfortunately to date no prediction is known for the spacing or the distribution of 
individual eigenvalues in the origin region relevant for chiral symmetry breaking. 
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5 Adjoint QCD and Symplectic Matrix Model 



In this section we summarise the known results for the (3d = 4 symmetry class of matrices with 
quaternion real elements, being relevant for SU(Nc) gauge theories in the adjoint representation. 
In the first subsection 15.11 explicit results for the complex eigenvalue correlation functions are given. 
Since the methods are similar to the previous section we will be brief here, and only highlight the 
main differences: the fact that for Nf ^ the densities remain real and positive as well as the different 
repulsion pattern of complex eigenvalues (see figs. [I] and [5]). The positivity makes both quenched and 
unquenched lattice simulations possible, and a comparison to our MM predictions is done in the next 
subsection 15.21 

5.1 Complex eigenvalue correlations from a symplectic MM 

Our starting point is the MM eq. (|4.16p with the matrices having quaternion real elements instead. 
The very same procedure described in the previous section can be employed to achieve a complex 
eigenvalue model. We refer to [51] for details of the derivation where this ensemble was first defined. 
The MM partition function for Nf flavours with equal chemical potential is reading in terms of complex 
eigenvalues 

J j=1 f=1 \ H- 

N N 

*U\4-#\*\4->? 2 \ 3 Il\4-tf\ 2 - (s-i) 

k>l h=l 

There are two main differences to eq. (|4.17p for [3d = 2. First, the mass term is automatically 
positive definite here because the eigenvalues of quaternion real matrices always come in complex 
conjugated pairs [H]. It thus resembles QCD with pairs of Nf conjugated flavours also called phase 
quenched (see fig. [7]). The same pairing of eigenvalues is also true for the zero modes, leading to a shift 
v — > 1v compared to QCD. Second, the Jacobian is different from the squared absolute value of the 
Vandermonde determinant for (3d = The last factor in eq. (15.11) leads to vanishing correlations 
along the x- and y-axis here, as \zi — ^ 2 | 2 = 16a;? y? vanishes there. We will find this feature explicitly 
in our results for the correlation functions, see fig. [TUl In the limit [i ~~ > for both (3d = 2,4 the 
Jacobians reduce to An(x 2 )^ d , counting the number of independent degrees of freedom per matrix 
element. 

The complex eigenvalue correlation functions are defined in complete analogy to eq. (|4.24|) with 
the corresponding weight and Jacobian. The result can be written as a Pfaffian |51[ [82] 

k 

R^( Zl ,...,z k ) = \[tz] 2 z])w( Zj ,z*) I'f,, , /, 

3=1 

For antisymmetric matrices A of even dimension the Pfaffian is the square root of the determinant, 
Pf(-A) = J det(j4). A similar form uses a quaternion determinant, see [82]. Here we have introduced 
the pre-kernel Kjv(zi, zf), 

N-l 1 

i^n(zi,4) = Y] — (q2k+i(zi)Q2k(4) ~ Q2k+i(4)l2k(zi)) , (5.3) 

— » ri. 



1 There exists a different sympiectic MM of normai matrices with Jacobian |A{2}| 4 , see e.g. [72] . 



KN{Zi,Zj, 
K N ( Z *,Zj 



K N ( Zi ,Z*) 
K N (Z*,Z*) 



(5.2) 
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containing skew orthogonal polynomials q k {z). These are defined by the following skew-orthogonality 
relation 

d 2 z w(z,z*)(z* 2 - z 2 ){q 2k+1 (z)q 2 i(z*) - q 2k+1 (z*)q 2l (z)} = r k S kl , (5.4) 



where the remaining integrals with even/even and odd/odd indices vanish. For the quenched weight 
with Nf = above these polynomials are again given in terms of Laguerre polynomials in the complex 
plane [51], 

-2 \ '■■ / \ - 2 , 

(5.5) 

j=0 



f ^ t2u ( Nz 2 \ * / Nz 2 \ 

q2k+l(Z) ~ L 2k+1 \^ 1 _ fx 2 j ' ~ 2^ c jW L 2j ^ _^ 2 J 



with some known /x-dependent coefficients Cj(fi) [51). The unquenched correlation functions can again 
be expressed in terms of the quenched pre-kernel |5H [82], where for an odd number of Nf flavours 
also the quenched skew orthogonal polynomials of even order q 2k (z) explicitly appear [82J. 

The large- iV limit is defined by the same rescaling of the chemical potential, eq. (|4.28p . of the 
complex eigenvalues and masse 

£3, ^2z = £ and ^2 m = rj, and correspondingly of the microscopic 
density eq. (|4.29p . The simplest examples for the resulting correlation functions in the large- N limit 
is the quenched microscopic spectral density 

pf U \0 = ^i(t 2 ~f) IC! 2 K 2u (g) e + ^« 2 ^ 2 )j,(C,D , (5-6) 

where the limiting pre-kernel eq. (|5.3f) leads to a double integral here due to the extra sum in eq. 
(f531) (compared to eq. (jOO]) for (3 D = 2) 



^(e,r)^^ 1 ^^ 1 ^e- 2s ( i+ ') a2 (j 2 ,(2 V / ^o^(2 v / ^r) - ^(2^0^(2^*)) . (5.7) 

Its antisymmetry together with the pre- factor (£* 2 — £ 2 ) leads to a real positive density eq. (|5.6p . 

Our second example is the unquenched microscopic spectral density for Nf = 2 flavours of degen- 
erate masses: 

Ps (0 " Ps (0 v u^m^M) ) ' (5 - 8) 



where 

K(Z,0 = J ds j dt J*j e ~ 2s ^ a2 (V~tJ 2l/+l (2V7t 0J2»(2V^ C) ~ J 2 „(2v^i t)<W(2^ 0) 

(5.9) 

is the derivative of the kernel eq. (|5.7p . It appears because of the degenerate masses in our example. 
Eq. (|5.8p now contains 3 terms being a Pfaffian, and it is again explicitly real. 

The quenched and unquenched densities eqs. (|5.6|) and (|5.8|) are compared in fig. [TU] for different 
values of rescaled chemical potential a and mass rj. The vanishing correlations along x- and y-axis due 
to the Jacobian is clearly observed. The introduction of masses has a mild effect, pushing the density 
further away from the origin when lowering the mass, without drastic change of its shape. At zero 
mass the resulting density becomes indistinguishable from the quenched density at v = 2. However, 
this flavour-topology duality holds only approximately for // ^ due to the explicit ^-dependence of 
the JT-Bessel weight. 



12 The reason for this difference is that both ensembles in eqs. (|5.ip and (|4.17p have the same variance at /i = 0, instead 
of ex.p[-/3 D x 2 /2]. 
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Figure 10: The quenched density eq. (|5.6p (top left) vs. the unquenched density eq. (|5.8j) for iVj = 2 
flavours of equal rescaled mass rj = 4.26 (top right), both at v = for a = 0.4, and the unquenched 
density eq. (15. 8p for Afy = 2 at a = 3.28 and rescaled mass 77 = 19 (bottom right). For comparison we 
also give the quenched density of real eigenvalues (bottom left) for v = (left curve) and v = 2 (right 
curve), respectively. 



For large values of a the density becomes almost constant along a strip, apart from the repulsion 
from the axis. In fig. [10] bottom right the mass at 7/ = 19 deforms this strip, as the location of the 
mass is a zero of the partition function, and of the density as well. The described effect of the masses 
on the density will be used in the next subsection to test the effect of unquenching in Lattice result. 

Finally let us give the microscopic density for an arbritrary number of flavours: 



w - ^ - ^ 2K - ) n if? + ' > <«■«» 

where the additional two flavours have masses i£ and The structure is exactly the same as for the 
unquenched densities for QCD, eq. (|4.38[) . The difference here is that due to the quaternion structure 
the extra complex mass pair £ and £* does not spoil the reality and positivity of the corresponding 
partition functions. For more details and higher order correlation functions we refer to \51\ 182] . 

Let us close this subsection with a few remarks. For (3d = 4 the complex eigenvalue correlation 
functions are simpler than the corresponding ones for \i = 0. This seemingly paradoxical statement 
originates in the simpler skew product and pre-kernel compared to the MM of real eigenvalues, ex- 
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plaining features of the latter [82] . 

In the previous section on the (3d = 2 class several alternative and equivalent derivations of the 
same results have been stated. In our (3d = 4 class to date no such alternative derivation is known 
so far using replicas, even for fj, = 0. Also the supersymmetric method has only been applied to 
non-chiral MMs with (3d = 4 (and (3d = 1) 02], where the microscopic density was derived. For a 
complete solution from skew-orthogonal polynomials of this non-chiral MM we refer to [85]. As in the 
(3d = 2 case a supersymmetric treatment would be simpler to do first in the one-MM eq. (|2,5p . 

Consequently little is known about the universality of the microscopic correlation functions and 
their relation to field theory in this symmetry class (the same statement is true for [i = 0, except for 
all masses degenerate). 

The only exception is for one pair of complex conjugate flavours Nf = 2, where the same result 
for the partition function 

^S 2 '%) ~ fdte-^I^tn) (5.11) 

can also be derived from the corresponding group integral in the exPT |88j. For more general MM 
partition functions with Nf > 2 we refer to [SU [82] . 




Reft) lm(9 

Figure 11: Lattice data (histograms) at mass ma = 0.07 in lattice unit a, volume V = 6 4 and fj, = 10 -3 
from [91 j . vs. MM prediction eq. (|5.8p at a = 0.0127 (full curve). The cuts are made parallel to the 
x-axis along the maxima (left), and perpendicular to it on the first pair of maxima. For comparison 
the quenched curve eq. (|5.6p is included as well, which is clearly shifted away from the data (left plot). 



5.2 Comparison to unquenched Lattice data using staggered fermions for SU{2) 

In this subsection we compare the complex eigenvalue correlation functions from the previous sub- 
section to lattice simulations [901 [9T] . Because for eigenvalues in the complex plane much statistics 
was needed staggered fermions were chosen in these simulations, using the code of [92J. In order to 
remain in the (3d = 4 symmetry class while using staggered fermions this comparison is made between 
our symplectic MM results and an SU(2) gauge theory in the fundamental representation realised 
as a staggered operator (see tabled]). We recall that in the continuum SU(2) in the fundamental 
representation belongs to the class (3d = 1- 

The lattice simulations [91] were done using two different volumes, V = 6 4 and 8 4 at fixed gauge 
coupling (3 = 1.3 for chemical potentials ranging from n = 10~ 6 — 0.2. Starting at large masses where 
the smallest eigenvalues effectively quench |9Q| the masses were varied down to values where an effect 
of unquenching was observed. For more details of the simulations we refer to [91J. 
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For a quantitative comparison we show cuts trough the 3 dimensional density plots of fig. |T0j the 
curves are cut along the maxima parallel to the x-axis and on the (first) maximum perpendicular to 
that parallel to the y-axis. 

For small a < 1 the effect of unquenching is only seen in the cut in x-direction through a shift of the 
curve, see fig. [11] left. The shift compared to the quenched curve from eq. (|5.6p is clearly visible. Note 
that it is very difficult to perform simulations at sufficiently small mass. In the perpendicular direction 
which is used to fit a no difference to the quenched curve is seen (when appropriately normalised) as 
we expect. The data follow very nicely the MM prediction. A similar comparison was made for large 
quark masses in [90] . nicely matching with the quenched result. 

Turning to large values of a the situation is reversed. At the smallest mass values available the 
cut in x-direction is not distinguishable from quenched. However, the perpendicular cut in y-direction 
shows the deformation of the density through the mass, as is clearly seen in fig. [T2l Since in this case 
the cut in y-direction is also used to fit a, two different fits are compared: the a determined from the 
unquenched curve is compared to the quenched curve at the same value of a, see the right curve in fig. 
[T2l Second, a different value of a. is obtained by directly fitting the quenched curve, leading however 
to a bigger mismatch in fig. [12] left, underscoring the data. Clearly the unquenched curve describes 
the data best. 
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Figure 12: The same cuts as in fig. [11] but for larger a: Lattice data (histograms) at mass ma = 0.06, 
volume V = 6 4 and \i = 0.2, from [91], versus MM prediction eq. (|5.8|) (full curve). For a better 
resolution we only display half of the cut in y-direction compared to fig. [TT] (right). Here two different 
quenched densities from eq. (|5.6p are inserted: the with the same a as unquenched, and a fit to the 
right cut. Both deviate from the data. 

The same procedure has been repeated for different values of the mass and chemical potential in 
[91j. Comparing the two volumes and keeping the product Vfj? fixed confirms the scaling predicted 
by e\FT and MMs from the data in the (3d = 4 symmetry class. 

In previous Lattice simulations using staggered fermions for both classes (3d = 4 and 1 |93j no 
comparison to MMs was available, and still is not for the latter class. 

In summary we have seen that our symplectic MM is able to describe unquenched Lattice data in 
the case where the density remains real and positive. 
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6 Matrix Model of QCD with Imaginary Isospin Chemical Potential 



In this section we discuss a MM for imaginary isospin chemical potential. The main difference to the 
previous MM is that is has real Dirac operator eigenvalues. This model was first introduced in [55] 
in terms of the corresponding e%PT. The motivation was to have a model that allows to determine 
the pion decay constant F n on the Lattice while keeping the Dirac eigenvalues real, in order to allow 
for unquenched simulations. The quenched results in [53] were then extended to the unquenched 
two-point function with Nf = 2 and confirmed on the Lattice [56J. We note that F n can also be 
determined on the lattice for real \x using quenched QCD with complex Dirac eigenvalues [57} I84j. 

Here we follow the presentation of [62] (see also |58| ) where the corresponding MM was introduced 
and solved for all unquenched correlations functions. The cases known from field theory [55} [56] were 
confirmed. 



6.1 Real eigenvalue correlations 

We start by writing down the Dirac operators for imaginary isospin chemical potential 

«"> = (-<W) *?•)■ (61) 

For both signs it is anti-hermitian thus having real eigenvalues. The symmetry classification in table 
Q] is not affected by this rotation and thus we can immediately write down the MM for m copies with 
-\-ijjL and n copies with — %[i 



/m n 
d®dm J] det[0 + + m f+ ] H det[0_ + m/_] e ~ N H^*+*^) . (6.2) 
-f i i -f 1 



/+=i f-=i 



However, for /j ^ the operators 0+(/i) and 0_(/x) now have different real eigenvalues (or more 
precisely real positive singular values): {x^} and {y^} G M + respectively. After the change of variables 

$± = $±/xtf (6.3) 

the Gaussian weight becomes replaced by 

V($+, $_) = c+($5.$+ + - c_ + (6.4) 



where c± = (1±^ 2 )/ (4^ 2 ). A transformation to eigenvalues of this two-MM can thus only be performed 
for the (3d = 2 symmetry class including QCD, to which we restrict ourselves in all the following. Only 
here the following matrix integral is available [M] 

dUdVexp[c-NTr(VXUY) + h.c. ] = — ^^f^A- — , (6.5) 

If we parametrise = U+XV+ and <£_ = U-YV- and use the invariance of the Haar measure we 
obtain the following eigenvalue model: 

poo N j m n 

Z^\{ mi },{m 2 }) ~ / [dx i dy i {x i y i y +1 (x 2 + m} + ) J] (y 2 + m}_) 

Jo i=i \ f+=i f-=i , 

xA J v({ a; 2 })A i v({y 2 })det[/ ! ,(2 C _iVx fc yO]^ Efc+{:!;?+ ^ ) . (6.6) 

k,l 



31 



While for the partition function the determinant can be replaced by its diagonal part due to the 
antisymmetry of the Vandermonde determinants, this is not the case for the eigenvalue correlations. 
We note that also in this model the initial Gaussian weight has become non-Gaussian containing an 
7-Bessel function, due to the group integral performed. 

The MM eq. (|6.6p can be completely solved by introducing bi-orthogonal polynomials 



dxdy wW(x,y) P^\x 2 ) Q^V) = 5 nk , (6.7) 
with respect to the weight function 

m n 

wW(x,y) = (xy) u+l 11 (x 2 + m 2 f+ ) J[ {y 2 + m)_) I u (2c^Nxy) e~ Nc+( - x2+y2) . (6.8) 

/+=i /-=l 

The monic polynomials are in general different. Together with their integral transform 

xl Nf \y) = rdxw( N f\x,y)P^\x 2 ), 



x[ Nf \x) = dxw^f\x,y)Q { k Nf \y 2 ). (6.9) 



four different kernels can be constructed 



K^\ y ,x) = y: Qk {y) h k — , h^\ x1 , X2 ) = y: Xk { i k ( a) 

k=0 k k=0 k 



i Nf \x,y) = E ~ I" ~ . H^\ yi ,y 2 ) = Y Qk {Vl)Xk iV2) 

k=0 k k=0 



hu — ' hh 

K k=o K 

These are needed to express all correlation functions of mixed x^ and yi eigenvalues defined as 
R 



? (N f ,v) N\ 



{n , k) ({*}„;{»}*) - (N-n)\(N-k)\4»f\ {mi}t{m2}) 

x / TT dxi U ^ det ™ {Nf \x P ,y q ) &N({x 2 })A N {{y 2 }) (6.11) 



i=n+l j=k+l 

H^ f ( Xil ,x i2 ) M { ^ f \x h ,y j2 ) - w( N f\ Xil ,y j2 ) 



det 

Kil,t 2 <n; l<7ij2<fc 



K {Nf) (v x- ) H {Nf) (v v 



(6.12) 



In the definition in the first line we have written the integrand of eq. (|6.6p in a more compact form. 
The solution in terms of a determinant of 4 blocks follows from [95J , adopted to our chiral two-MM 
|62j . Note that in contrast to eqs. (|4.25p and (|5.2p the kernels built from the integral transforms of 
the polynomials contain part of the weight inside the determinant. 

For Nf all polynomials, their transforms and the corresponding kernels can be expressed 
through the quenched quantities as detailed in [62]. Below we give the quenched quantities as an 
explicit simple example. 

For the quenched weight Nf = in eq. (|6.8p the bi-orthogonal polynomials are equal and both 
given by Laguerre polynomials of real variables: 

P n (x 2 ) , Q n (x 2 ) ~ L v n (-^xA • (6-13) 
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In addition the integral transforms also become proportional to Laguerre polynomials, times the 
Laguerre weight: 

Xk(y) , My) ~ y 2v+1 e ^ y L\ W^V 2 ) • (6-14) 

The reason is that for this special case the bi-orthogonal polynomials and their weight can be con- 
structed starting from two ordinary orthogonal Laguerre polynomials of a factorised weight, where for 
details we refer to appendix B of ref. [62J . 

We finish this subsection by giving two examples in the microscopic large- N limit given by eq. 
(|4.28|) . First, the partition function eq. (|6.6|) in this limit becomes identical to eq. (|'2.12p . with the 
replacement fi 2 — > —fi 2 in eq. (|2.13p [62J. It agrees with the partition functions computed from e%PT 
field theory, simply because the effective chiral Lagrangian eq. (12. 8p is obtained by this rotation to 
B = ifj, diag(l| m , — l| n ). The group integral is convergent for both signs of fi 2 and the rotation is 
trivial. 

On the level of correlation functions however, this rotation is highly nontrivial. Leading to complex 
eigenvalues on the one side and two sets of real eigenvalues on the other side their determinantal 
structure is very different, comparing eqs. (|4.25p and (|6.12p above. 

The simplest nontrivial correlation function is the quenched probability to find one eigenvalues £ 
of Ip+ and one eigenvalues £ of J])—. In the microscopic limit eq. (12. 9p we obtain from eq. (|6.12p 

plT$(fcO - pf f ' U \0 pP"\0 + pf^^O ■ (6.15) 
Here the first product contains the quenched microscopic density eq. (I4.3ip at zero 

pf^ho = |[J2(o-jH-i(OJi-i(0] , 

whereas the connected part is given by 

Here the same building block as in eq. (|2.13p appears (see also eq. (|4.30p ). containing the real 
eigenvalues 

T±(£,C) = f dtte ±2 » 2 *UZt)MCt) ■ (6.17) 
J o 

The connected density is plotted in the next subsection, fig. CGI 
6.2 Comparison to Lattice data 

Because the introduction of an imaginary isospin chemical potentials does not change the anti- 
Hermiticity of the Dirac operator, standard Monte-Carlo simulations are possible in the unquenched 
case as well. In this section we give a very brief account of the lattice results of [55] [56] to where we 
refer for details. There, staggered unimproved fermions were used, which does not lead out of the 
/3d = 2 symmetry class. The simulations covered several volumes V = 6 4 , 8 4 and 10 4 at various small 
values of isospin chemical potentials [ii so = 0.002 — 0.01. 

In order to better understand the pictures below let us rewrite the correlation function eq. (|6.15p 
used there: 

pf$ ) (t = 2Nx;t = 2Ny) = Vim^ <V S\x - x k ) £ 8\y - Vj ) \ (6.18) 
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The connected part that will be compared to data is obtained by subtracting the factorised correlation 
function in eq. (|6. 15j) . For /i ^ the eigenvalues of the two operators lft±(n) are clearly distinct. In 
the limit fi — > however, they will coincide, leading to a delta-function in the connected correlator 
~ 5 1 (y — x). The effect [i 7^ has is to smooth out this delta- function, with the width of the resulting 
curve being sensitive to F n contained in the rescaling factor of /i. In fig. [13]left the quenched connected 
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Figure 13: Quenched (left) and unquenched [SB] (right) Lattice data at V = 8 4 (left) and a = 0.152 
at V = 6 4 (right), vs. the connected correlation Pgf^i\ C ° nn (£,C) (see eq. (]6. 16|) for quenched). 

P^(i'l) COnn (£'0 density eq. (16. 16ft is plotted versus lattice data [55]: here the second argument is 
kept fixed at £ = 4 while the first argument £ is varied. For comparison the same correlation function 
in the limit fi — > is inserted in dashed line, omitting the delta-spike at £ = £ = 4. 

In the second plot fig. [13] right the unquenched density generalising eq. (|6.16j) (see [56] [62] for 
details) is compared to the data. The data nicely fit the curve corresponding to bare mass m = 0.002 
(lowest curve to the left), distinguishing well from the parameter set inserted for comparison m = 0.005 
(see [56] for the corresponding simulations). The quenched limit obtained by setting 2Nm = oo is 
inserted for comparison as well (lowest curve on the right). This example underlines once more the 
fruitful interplay between MMs, effective field theory and Lattice simulations. 



7 Conclusions and Open Problems 

We have seen that Matrix Models are a useful tool to investigate QCD with both real and imagi- 
nary chemical potential, extending their previous success for zero T and zero \x. MMs provide both 
qualitative predictions such as for the phase diagram, and quantitative results for the Dirac operator 
spectrum that can be confronted with Lattice simulations. 

Following the symmetry classification into 3 classes we have explicitly shown how to construct and 
solve MMs with complex eigenvalues for two of these 3 classes, the class with three or more colours 
Nc > 3 containing QCD, and the class of adjoint QCD. All complex eigenvalue correlation functions 
could be computed for these two classes as a function of the number of quark flavours Nj, their masses 
wif and chemical potential of equal modulus ±/i. In particular this included unquenched correlation 
functions in QCD where the spectral densities become complex valued quantities. 

The same construction was then repeated to build a MM for the QCD class with imaginary isospin 
chemical potential. Here all real eigenvalues correlation functions were computed for the two sets of 
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Dirac operators lfi(±i[i). 

Several issues have been understood on the way. The equivalence with the underlying effective 
field theory, the epsilon regime of chiral perturbation theory was shown for the QCD class both for 
real and imaginary /x, on the level of the partition function for an arbitrary set of chemical potentials 
Hf. We have seen for QCD why the partition function with ordinary quarks and the corresponding 
sum rules remain ^-independent, while the generating functional for the density is //-dependent. The 
strong oscillations and the complex nature of the unquenched eigenvalue density in QCD have been 
made responsible for chiral symmetry breaking, in contrast to quenched and phase quenched QCD. 

The detailed MM predictions for the microscopic density explained the different signatures of 
complex Dirac eigenvalues, including their attraction or repulsion from the real and imaginary axis. 
Excellent agreement was found when comparing these quantitatively to Lattice simulations. We made 
comparisons to quenched QCD using staggered and overlap fermions, allowing to compare to different 
topological sectors of the gauge fields. For the adjoint QCD class we could compare to unquenched 
simulations. Here staggered fermions with two colours and fermions in the fundamental representation 
were used, which are in the same symmetry class. Again excellent agreement was found when lowering 
the mass to effectively unquenched the eigenvalues close to the origin. 

For imaginary isospin chemical potential in the QCD class we could make the same comparison 
between analytic predictions for two Dirac operators with different real eigenvalues and Lattice data 
using staggered fermions. Both for quenching and unquenching an excellent agreement was found. 

Several interesting problems remain open within the MM approach. Fortunately the symme- 
try class including unquenched QCD is technically the simplest within MMs. Nevertheless partially 
quenched correlation functions with several different \ij have so far only been computed for imaginary 
chemical potential for 2 flavours, fii / ^2- It would be very interesting to extend this to more flavours, 
and in particular to real fi, allowing for instance to compute correlations with in a background 

of zero chemical potential. 

To match with the partially quenched we would have to compute the corresponding group 
integral for a general baryon charge matrix B. Although we have shown that the expression for MMs 
and xPT are equal, we have not yet been able to compute the integral in either way. 

Turning to the 2 other symmetry classes only the class for adjoint QCD corresponding to symplectic 
MMs has bee completely solved. Here the link to e^PT is an open question for Nf > 2. The 
(3d = 1 class for two-colour QCD remains the most challenging problem for MMs. Even the easier 
corresponding non-chiral model of real non-symmetric matrices has resisted a solution for 40 years. 
This class would be very interesting also from a Lattice perspective as it has a sign problem: although 
the determinant det[j0(n) + m] remains real its sign fluctuates. 

Lattice simulations of these 2 latter classes may now be feasible with chiral fermions at non- 
vanishing fi due to the existence of the corresponding overlap Dirac operator at (3d = 2, without 
having to use the staggered formulation. Such simulations would be very interesting to perform. The 
most exciting challenge of course remains to see the MM for unquenched QCD realised on the Lattice. 

As a final remark it is well known that a much larger class of MMs with complex eigenvalues exists, 
compared to the 3 chiral classes for fi = 0. It would be very interesting to identify and solve some of 
the other classes that may be relevant in wider applications as well. 
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